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ABSTRACT 


The hot gas that constitutes the intracluster medium (ICM) has been studied at X-ray and 
millimeter/sub-millimeter wavelengths (Sunyaev-Zeld'ovich effect) for decades. Fast radio bursts 
(FRBs) offer an additional method of directly measuring the ICM and gas surrounding clusters, via 
observables such as dispersion measure (DM) and Faraday rotation measure (RM). We report the 
discovery of two FRB sources detected with the Deep Synoptic Array (DSA-110) whose host galaxies 
belong to massive galaxy clusters. In both cases, the FRBs exhibit excess extragalactic DM, some of 
which likely originates in the ICM of their respective clusters. FRB 202209144 resides in the galaxy 
cluster Abell2310 at z = 0.1125 with a projected offset from the cluster center of 520450 kpe. The 
host of a second source, FRB 20220509G, is an elliptical galaxy at z = 0.0894 that belongs to the 
galaxy cluster Abell 2311 at projected offset of 8704-50 kpc. These sources represent the first time an 
FRB has been localized to a galaxy cluster. We combine our FRB data with archival X-ray, SZ, and 
optical observations of these clusters in order to infer properties of the ICM, including a measurement 
of gas temperature from DM and ysz of 0.8-3.9 keV. We then compare our results to massive cluster 
halos from the IllustrisTNG simulation. Finally, we describe how large samples of localized FRBs from 


future surveys will constrain the ICM, particularly beyond the virial radius of clusters. 


Keywords: galaxy clusters — radio transient sources, Abell 2310, Abell 2311 


1. INTRODUCTION 


Galaxy clusters are massive (10'4 — 10? Ma) gravita- 
tionally bound objects comprised of hundreds to thou- 
sands of galaxies. Galaxies make up just a few percent 
of the total cluster mass. The dominant component of 
baryons by mass and volume is in the hot intracluster 
medium (ICM) made of diffuse gas with n, ~ 107? cm? 
and T, ~ 10" — 108 K. Within a sufficiently large vol- 
ume (e.g., the virial radius), the ratio of dark matter 
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to baryons approximately matches the Universal value 
(Eckert et al. 2019), which is not true for smaller ha- 
los where feedback is expected to expel gas beyond the 
virialized dark matter halos (Tumlinson et al. 2017). 
The ICM has been studied in great detail at X- 
ray wavelengths for the past fifty years (Sarazin 1986). 
Thermal bremsstrahlung emission from gas heated to 2— 
15 keV is observed out to the virial radii of clusters (1— 
2 Mpc), tracing the radial distribution of the ICM gas 
(Ghirardini et al. 2019). The ICM reaches such high 
temperatures by adiabatic compression and shock heat- 
ing as the gas reaches hydrostatic equilibrium within 
the potential well (Voit 2005). Free-free emission scales 
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with emissivity, and the specific X-ray luminosity can 
be written as an integral of plasma density squared, 


Ly = | n &(r)av. (1) 


where ne is the number density of free electrons and 
A(T.) characterizes the temperature dependence of X- 
ray emission. The total mass can then be inferred from 
the cluster's X-ray luminosity (Pratt et al. 2009a). X- 
ray spectroscopy can also be used to measure the gas 
temperature (Bóhringer & Werner 2010) and its velocity 
structure (Hitomi Collaboration et al. 2016). 

The hot plasma in the ICM is also observable 
at millimeter/sub-millimeter wavelengths through the 
thermal Sunyaev-Zel'dovich (SZ) effect, whereby pho- 
tons in the cosmic microwave background (CMB) 
are inverse Compton scattered by thermal electrons 
(Mroczkowski et al. 2019). This leads to spectral distor- 
tions in the CMB at the level of 1074-1075. SZ obser- 
vations of the ICM are parameterized by the Compton 
Y-parameter, which is given by the integral of electron 
thermal pressure along the line of sight, 


k 
Uysz = Bos [n.i (2) 
Me C 


where kg is Boltzmann's constant, c is the Thomson 
cross-section, Me is electron mass, c is the speed of light. 

Fast radio bursts (FRBs) offer an additional probe of 
the ICM of massive galaxy clusters. FRBs are short- 
duration (107?—10? ms) radio transients that have been 
detected over cosmological distances (z S 1.5) (Petroff 
et al. 2019; Cordes & Chatterjee 2019). To date, O(10?) 
sources have been discovered (CHIME/FRB Collabora- 
tion et al. 2021; Petroff et al. 2016) but only two-dozen 
have been localized with sufficient angular precision to 
identify a host galaxy (Chatterjee et al. 2017; Ravi et al. 
2019; Heintz et al. 2020; Bhandari et al. 2020). Propaga- 
tion effects imparted on the radio pulse encode informa- 
tion about the plasma through which the burst travelled. 
Dispersion measure (DM) is an integral of line-of-sight 
electron density and is given by, 


DM = I nodi: (3) 


Unlike with thermal X-ray emission or SZ observations 
of the ICM, the FRB DM is detectably impacted by all 
plasma between us and the source, not just 10975 K gas. 
'This is both a blessing and a curse: FRB DMs will probe 
plasma in the host galaxy, the circumgalactic medium 
(CGM), the ICM of intervening or host clusters, the 
intergalactic medium (IGM), and the Milky Way's in- 
terstellar medium (ISM). For individual FRB sightlines, 


these terms must be explicitly modelled. An FRB that 
is impacted by a galaxy cluster will have the following 
“DM budget”, 


DMicm , DMhost 
1+ 2 | 1+ Zfrb ' 


DMops = DMyw + DMiem 4 


where z, and zr, are the cluster and FRB redshifts, re- 
spectively. The observed DM is likely dominated by the 
IGM for most FRBs. Still, it is important to be specific 
about what is meant by the IGM. In Macquart et al. 
(2020), the relationship between DM osmic and source 
redshift includes all extragalactic gas between the Milky 
Way and the host galaxy. DMcosmic therefore includes 
the CGM, intragroup medium (IGrM), ICM, and IGM. 
All of those terms correlate with source distance because 
the optical depth of halos increases with redshift. In this 
work we take the IGM to be ionized gas that resides out- 
side of virialized dark matter halos (McQuinn 2016). 

For sources that intersect or are embedded in galaxy 
clusters, the largest DM contribution may be from the 
ICM (Prochaska & Zheng 2019). For example, a cluster 
with (ne) = 107? cm^? and Roo = 1 Mpc will lead to 
DMicm ~ 10? pe em"? for low impact parameters. This 
value is more than double the typical contribution from 
the IGM for a source at z — 0.5 and 7096 larger than the 
mean DMoss of CHIME/FRB sources (Macquart et al. 
2020; CHIME/FRB Collaboration et al. 2021). The po- 
larization properties of FRBs can help measure the mag- 
netic field of the ICM. Faraday rotation measure (RM) 
is an integral of the free electrons, weighted by the line- 
of-sight magnetic field strength, 


3 
RM = 


2T m2 ct 


where e is the charge of an electron. While other ra- 
dio objects can provide cluster RMs (van Weeren et al. 
2019), FRBs are unique amongst extragalactic sources 
in their ability to measure both RM and DM. Therefore, 
if one can isolate the components of RM and DM that 
are due to the ICM, the mean line-of-sight magnetic field 
strength in the cluster can be determined by their ratio, 


RMicm (6) 
DMicm’ 


where RM and DM are given in the standard units of 
rad m^? and pecm~°, respectively. The application of 
FRBs to the ICM has been considered in previous works, 
including the detection of an FRB in the direction of the 
Virgo cluster by ASKAP (Agarwal et al. 2019). In that 
case, the localization precision was insufficient to deter- 
mine a host galaxy, so it was unclear if the source was 
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embedded in, or behind, the cluster. Another unlocal- 
ized source, FRB 20190116A (CHIME/FRB Collabora- 
tion et al. 2019) is likely within a ~ 5 Mpc of the Coma 
cluster in projection and may be dispersed by filamen- 
tary structure extending from the cluster (Hallinan et al. 
2019a). Fujita et al. (2017) investigated the application 
of FRBs to the warm hot intergalactic medium (WHIM) 
at the outskirts of galaxy clusters, where X-ray obser- 
vations are difficult due to the n2 dependence of emis- 
sivity. Nonetheless, most of the effort in studying halo 
gas with FRBs has been devoted to the circumgalac- 
tic medium (Prochaska & Zheng 2019; Ravi 2019; Lee 
et al. 2022). That is because the CGM is more difficult 
to detect than the ICM using traditional means (X-ray 
and SZ observations are much less sensitive to the gas 
in galaxy-scale halos), and because of the CGM’s signifi- 
cance in galaxy formation and feedback processes. Con- 
nor & Ravi (2022) presented the first statistical evidence 
for the impact of halo gas on FRBs and demonstrated 
the importance of considering galaxy groups on FRB 
DM budgets. As we demonstrate, current and future 
FRB surveys will find a significant number of sources 
that are impacted by galaxy clusters, and understand- 
ing their effects and prevalence are important for CGM 
studies, in addition to studying the ICM itself. 

In this work we make the first unambiguous detection 
of FRB sources that belong to a galaxy cluster. Both 
sources were discovered by the real-time FRB survey on 
the Deep Synoptic Array (DSA-110)!. The DSA-110 is 
a radio interferometer at the Owens Valley Radio Ob- 
servatory (OVRO) that was built to localize FRBs with 
~ arcsecond precision. A detailed account of the instru- 
ment will be presented in Ravi et al. (in prep.). During 
the observations presented in this work, the DSA-110 
was operating in its commissioning phase with 63 anten- 
nas collecting data, as described in Ravi et al. (2022). 

In Section 2 we describe the radio properties of two 
FRBs using the DSA-110. We describe the two massive 
clusters of which the FRB host galaxies are members. 
A companion paper details the interferometric localiza- 
tions of FRB 20220509G and FRB 20220914A as well as 
their unusual host galaxies (Sharma et al, in prep.). 
In Section 3.2 we model the DM contribution from the 
ICM for both sources, combining our radio data with 
archival X-ray, SZ, and optical data to infer properties 
about the clusters’ respective ICM. Finally, we compare 
our results with clusters found in the IllustrisTNG sim- 
ulation. 


1 https://deepsynoptic.org 


2. CLUSTER FRBS 


The radio properties of FRB20220914A and 
FRB 20220509G are listed in Table 1. The observed 
parameters of their respective host galaxy clusters are 
given in Table 2. The two FRBs presented here are 
from distinct galaxy clusters, but are relatively nearby 
on the sky—separated by just 3.1 degrees. DSA-110 
is a transit instrument that was parked at Declination 
~ +70 for most of the 2022 commissioning period. For 
this reason, most of our initial sample of sources fall on 
a ring in RA and several pairs are within a few degrees 
of each other. Below we described observations of the 
FRBs and their host galaxies, as well as archival data 
obtained for the galaxy clusters. 


S/N 


Frequency (MHz) 


-0.5 0.0 0.5 
Time (ms) 


S/N 


Frequency (MHz) 


-15 -10 -0. 0.0 0.5 
Time (ms) 


Figure 1. Dedispersed total-intensity dynamic spectra for 
two FRBs detected by DSA-110. Both FRBs reside in galax- 
ies that belong to massive galaxy clusters. 


2.1. FRB 20220914A 


FRB 20220914A was detected at MJD 59836.1459660. 
Follow-up analysis found an optimal DM of 
631.3 pccm™? and a pulse full-width at half maximum 
(FWHM) of 140 us. With a S/N of 13.4 and system- 
equivalent flux density (SEFD) of 140 Jy, the fluence of 
the burst was roughly 2.2 Jy ms. Offline interferometric 
localization found the position of the FRB to be R.A. 
J2000, decl. J2000 = 18h48m13.63s, +73d20m12.89s, 
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Source RA (deg) Dec (deg) Zgal DM; DMuw RMəs T (ms) F (Jy ms) b, (kpc) 
FRB 20220914A 282.0568 73.3369 0.1139 631.29 47/55 N/A < 0.080 2.2 520+50 
FRB 20220509G 282.6700 70.2438 0.0894 269.53 46/55 —110(1) 0.08+0.02 5.5 870+50 


Table 1. The radio properties of two galaxy cluster FRBs. DM and RM are reported in pccm~? and rad m7’, respectively. 
The scattering timescale, 7, is referenced to 1.4GHz. The two Galactic ISM DM estimates listed are YMW16 and NE2001. RA 


and Dec are in epoch J2000. 


Cluster RA (deg) Dec (deg) Richness Zclust Ms00 (Mo) Rsoo (kpc) Lx (erg/s) Ysrso0 (arcmin?) 
Abell 2310 281.81708 73.3472 124 0.1125 1.7x10"* 810 8.3x 10 0.85 + 0.35 x 10? 
Abell2311 282.497 70.376 190 0.0899(1) — 1.6x10'* 800 4.7x10*% 0.81 0.36 x 10? 


Table 2. The properties of two galaxy clusters that host the FRBs presented in this paper. Lx is the X-ray luminosity at 


0.1-2.4 keV within Rsoo. 


the centroid of which is 1.6" from a galaxy with spec- 
troscopic redshift z — 0.1139. The expected DM from 
the Milky Way at these Galactic coordinates is roughly 
50pccm ? from NE2001 (Cordes & Lazio 2002) and 
47 pccm ? from YMWI106 (Yao et al. 2017). The IGM 
is expected to contribute 50-120 pccm-? at this red- 
shift, suggesting that FRB 20220914A has significant 
extragalactic DM excess. 

We find no evidence of scattering, placing an upper- 
limit of 7 < 80ys at 1.4GHz. Polarimetric analysis 
found no detectable polarization and no RM was deter- 
mined after searching a range of —10° to +10° rad m~?. 
The upper-limits on polarization fraction were L/I « 
1596 and V/I « 20%. 


2.1.1. ABELL 2310 galaxy cluster 


The position of FRB 20220914A was cross-matched 
with both the DESI Legacy Imaging Surveys Data Re- 
lease 9 (DR9) galaxy catalog (Dey et al. 2019) and DR9 
cluster/group catalog (Yang et al. 2021). We find that 
the host galaxy of FRB 20220914A is a member of the 
massive galaxy cluster Abell2310 (Abell 1958). Ac- 
cording to the DESI DR9 group/cluster catalog, this 
cluster has richness 124 (i.e. number of member galax- 
ies), Miso = 2.5 x 1014 Mo, with a brightest cluster 
galaxy (BCG) at spectroscopic redshift z — 0.1125. 
The cluster is also in the Meta-Catalog of X-Ray De- 
tected Clusters of Galaxies (MCXC) with source name 
J1847.2-- 7320 (Piffaretti et al. 2011). The X-ray surface 
brightness centroid is at 18h47m16.s 73d20m50s (Pif- 
faretti et al. 2011). The FRB offset from this position 
is 4.22’ suggesting a projected physical impact param- 
eter of 520kpc. The X-ray luminosity within R509 is 
8.31x10* ergs-!, which gives Msoo = 1.69 x 10!^ Mo 
using an empirical Lx — M relation (Pratt et al. 2009b; 
Piffaretti et al. 2011). 


'The cluster has also been detected via the thermal 
SZ effect and is listed in the Planck SZ2 cluster cata- 
log as PSZ2 G104.29--26.17 (Planck Collaboration et al. 
2016a). The SZ-derived mass based on a hydrostatic 
mass calibration is Msg9 = 1.94 + 0.28 x 10!^ Mo 
(Planck Collaboration et al. 2014), and its SZ centroid 
has a physical offset from FRB 20220914A of ~ 380 kpc. 
We have also analyzed the Planck MILCA y-map from 
Planck Collaboration et al. (2016b) using the techniques 
described in Sayers et al. (2016) to obtain both a two- 
dimensional projected model of the cluster y(0), along 
with the total integrated SZ signal Ys p59. In brief, we 
assume the cluster follows the profile shape given by Ar- 
naud et al. (2010), with the scale radius set by the X-ray 
value of Rso9 and centered on the X-ray centroid. We 
then determine the best-fit normalization of this model, 
with the uncertainty on this normalization estimated 
from fits to 100 random realizations of the Planck noise. 
The result is Y5nsoo = 0.85 + 0.35 arcmin?. 

The point-spread function (PSF) of Planck is large, 
so we take the offset between FRB 20220914A and 
Abell 2310 to be the separation between the X-ray cen- 
troid and the FRB's host galaxy. This gives a projected 
physical offset of 520-E50 kpc. In Figure 4 we show the 
filtered X-ray data with a contour for the SZ emission 
region, as well as a zoom in on the optical host galaxy 
image. The X-ray intensity data is taken from the public 
archive of X-ray clusters? (Xu et al. 2022). 

'The redshift difference between the host galaxy of 
FRB 20220914A and the brightest cluster galaxy (BCG) 
is Az — —0.00129, indicating the host has a recession 
velocity of roughly 390 km s^!. We do not know the ra- 
dial position of the FRB within the cluster, but we have 
good reason to believe it is not near the front. We take 


2 https:/ /github.com/wwxu/rxgcc.github.io 
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FRB 20220509 
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Figure 2. The polarization properties of FRB 20220509G. 
The top panel shows a flat PA across the pulse and the mid- 
dle panel shows the linearly polarized pulse profile. The 
bottom panel is the Faraday spectrum with a peak at 
RM--111.5rad m ?. 


the radial position to be lros, which is zero at the dis- 
tance of the cluster center and negative behind the clus- 
ter. T'he radial distribution of galaxies in clusters is well 
described by an NFW profile with concentration param- 
eter c=2.6 (Budzynski et al. 2012). With a virial radius 
of 800 kpc, the probability distribution of galaxy radius 
peaks around 300 kpc and declines towards larger radii. 
Therefore, for an impact parameter beyond 300 kpc, the 
most likely lzos is zero. More convincingly, as we show 
in Section 3.2, the FRB has significant excess disper- 
sion that likely comes from the ICM. The magnitude of 
this extra DM cannot easily be explained by the IGM 
or the host galaxy, so [ros is likely smaller than several 
hundred kpc. 


2.1.2. The host galazy of FRB 20220914A 


The FRB 20220914 has been localized to a typical 
late-type spiral galaxy at R.A. (J2000) = 18:48:13.9580 
and declination (J2000) = +73:20:10.703 (Sharma et 
al., in prep.). The spectroscopic redshift of the galaxy, 
as measured from optical spectra acquired with Keck- 
I/LRIS, is 0.1139 + 0.0001. A detailed spectral en- 
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Figure 3. Slice plots of the smoothed galaxy number density 
within 5 arcminutes of each FRB line-of-sight. The galaxy 
data are from the DESI Legacy Imaging Survey DR9 and are 
mostly photometric redshifts. 


ergy distribution analysis revealed a stellar mass of 
log M,(Mo) = 9.09* 503. 'The constrained star forma- 
tion history indicates significant recent star formation 
with several starbursts over the last 3.5 Gyr, thus re- 
flecting a wide probability distribution for the age of its 
progenitor. 


2.2. FRB 20220509G 


The fast radio burst source FRB 20220509G had an 
arrival time of 59708.4944991 at reference frequency of 
1500 MHz. Its optimal DM was 269.53 pccm™? and 
had a fluence of 5.5Jy ms. Offline interferometric lo- 
calization found that the source is at R.A. J2000, decl. 
J2000 = 18h50m40.8s, --70d14m37.8, with a 90% er- 
ror ellipse with axes 4.7 and 3.2 in R.A. and decl. 
This is within 6” of a massive galaxy with spectro- 
scopic redshift z = 0.0894. The Galactic ISM DM 
estimates are 55 pccm™? from NE2001 and 46 pc cm~? 
from YMW16. The expected IGM contribution at this 
redshift is roughly 75 pc cm-?, but with significant un- 
certainty. We analyze the source's DM budget in Sec- 
tion 3.2. The burst was fit with a scattering tail plus 
a Gaussian component. We find evidence of scattering 
with a timescale of 8020 us at 1.4 GHz and a Gaussian 
component that is 30 + 10 us. Offline polarimetric anal- 
ysis found that FRB 20220509G was nearly 100% lin- 
early polarized with an RM of —111.5 + I.5rad m^? in 
the observer frame. This is significantly larger than the 
Galactic RM foreground value of —9 + 16radm~*. The 
burst’s position angle (PA) is flat across the pulse. The 
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polarization properties of FRB 20220509G are shown in 
Figure 2. 


2.2.1. Abell2311 galaxy cluster 


The host galaxy of FRB 20220509G is a member of 
the galaxy cluster Abell2311 (Abell 1958). According 
to the DESI DR9 cluster/group catalog, the cluster's 
richness is 190 with Miso = 2.5 x 10'4 Mo, determined 
by the velocity dispersion of member galaxies. The ra- 
dius at which the average enclosed density is 500 times 
the critical density is R599 = 800kpc. The X-ray lu- 
minosity within R599 is Lx = 4.7 x 10“ ergs! with 
an inferred mass of Msoo = 1.6 x 10!^ Mc. The cluster 
does not have a published SZ detection and is not in the 
Planck 8Z2 cluster catalog. However, we have analyzed 
the Planck MILCA y-map for this cluster in the same 
way as for Abell 2311, finding Y5nsoo = 0.8140.36 10^? 
arcmin?. 

The fillered ROSAT X-ray intensity is shown in 
Figure 4. The projected physical offset between 
FRB 20220509G and the X-ray centroid is 870 + 50 kpc, 
placing the FRB at a minimum radius of just beyond 
Fisoo- 


2.2.2. The host galazy of FRB 20220509G 


The likely host galaxy of FRB 20220509 is an early- 
type elliptical galaxy with insignificant ongoing star for- 
mation, and hence stands out as a quiescent galaxy in 
the known population of FRB hosts (Sharma et al., in 
prep.). The measured spectroscopic redshift of the host 
is 0.0894 + 0.0001 and a detailed spectral energy distri- 
bution analysis reveals a stellar mass of log M. (Mo) = 
bl with star formation rate averaged over the 
last 100 Myr of 0.0870 ta: An old stellar population im- 
plies a long delay between the time of occurrence of this 
FRB and the formation of its progenitor, thus opening 
another window of progenitors with long delay times. 


3. DATA ANALYSIS 


We now seek to synthesize the FRB radio properties 
with observations of their host galaxies and host clus- 
ters. This will allow us to infer the impact of the ICM 
on FRB 20220509G and FRB 20220914A. Neither host 
galaxy has a high rate of star formation. The early-type 
host of FRB 20220509G, in particular, is not expected 
to have a significant ISM that would contribute to lo- 
cal dispersion (Lees 1992). In this section we estimate 
the origin of scattering, DM, and RM and model the 
contribution of the ICM to each propagation effect. 


3.1. Scattering & scintillation 


Both FRB sources show evidence of scintillation that 
is consistent with a nearby scattering screen in the 


Milky Way. FRB 20220509G and FRB 20220914A have 
decorrelation bandwidths of Av = 3.1 + 0.7 MHz and 
Av = 2.51 + 0.16MHz, respectively. Since the two 
sources are only separated by a few degrees, it is un- 
surprising that their decorrelation bandwidth and mod- 
ulation indexes are similar. However, the presence of 
fully-modulated scintillation is physically interesting in 
the case of FRB 20220509G, which also shows evidence 
of temporal scattering. The presence of both scattering 
and spectral scintillation allows one to place constraints 
on the geometry of the two scattering screens. Scatter- 
ing causes angular broadening and if that broadening 
were too large the FRB would be resolved out by the 
scattering screen in our Galaxy (Masui et al. 2015; Con- 
nor et al. 2016; Simard & Ravi 2021). 

For a given scattering timescale, angular broadening 
is maximal when the screen is halfway between us and 
the source. Therefore, if the screen were in the CGM 
of an intervening galaxy at z ~ 0.05, we would not ex- 
pect FRB 20220509G to scintillate. Following Eq 48 in 
Simard & Ravi (2021), we place an upper-limit on the 
distance between the FRB source and the first scattering 
screen to be < 140 kpc. This leads to three scenarios for 
the origin of the temporal scattering in FRB 20220509G, 
each of which is plausible but somewhat surprising given 
the early-type host galaxy and the hot, smooth ICM in 
which it is embedded. The scattering could be in the 
ISM of the host galaxy at ~ kpc scales from the FRB 
emitting source. This would require an unusual sightline 
given the lack of Ho emission and limited turbulence in 
the ISM of elliptical galaxies (Seta et al. 2021; Ocker 
et al. 2022). Alternatively, it could arise in the imme- 
diate vicinity of the source, such as in a stellar wind, 
analogous to the scattering in FRB 20190520 (Niu et al. 
2022a; Ocker et al. 2023; Beniamini et al. 2022; Anna- 
Thomas et al. 2022). A third option is that it could be 
near the host galaxy in the ICM, perhaps in the outflows 
or CGM of the host. 


3.2. DM from the ICM 


We model the DM along the line of sight to each 
FRB in order to generate a probability distribution for 
the ICM contribution. The probability density function 
(PDF, denoted here by P) of a sum of independent vari- 
ables is the convolution of their individual PDFs. The 
DM terms are not strictly independent, but we find that 
the convolution relation for a sum of variables is a good 
approximation in this case. Using Eq. 4, we find, 


P (DMicm) = P (DMobs)*P (-DMwuw)*P(-DMiau)*P(-DMi) 


(7) 
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Figure 4. The two FRB sources reported in this work with their respective host galaxies and galaxies clusters. The left column 
images show filtered ROSAT X-ray intensity (Xu et al. 2022) for galaxy cluster Abell 2310 (top) and Abell 2311 (bottom). The 
magenta circle in the top left panel has a radius of 1 Mpc and denotes the location of the SZ cluster from the PSZ2 catalog. 
The bottom panel magenta cirlce corresponds to our SZ measurement for Abell 2311. The right columns show insets on the 
position of FRB 20220914A (top) and FRB 20220509G (bottom) plotted over a PanSTARRSI riz-band image (Chambers et al. 


2016, PS1). 


where P(DMops) is Gaussian with variance determined 
by the measurement error on the observed DM. Here 
the DMs are all in the observer frame, so we must mul- 
tiply by 1+ ze when estimating the DM from the ICM. 
We assert that the probability density must be zero for 
negative DM. 


For FRB 20220914A, we place an upper-limit on DM), 
to be 50pccm-? based on the strong upper limit on 
temporal scattering. We take P(DMy) to be a uni- 
form distribution between 0 and 50 pc cm"? because we 
do not have a strong motivation for a preferred DM 
peak in the distribution. Previous detections of un- 
scattered FRBs do not have a robust measurement of 
local DM on which to base an empirical distribution. 
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There is reason to believe the local DM contribution 
is very small for other similar FRB sources. For ex- 
ample, FRB 20180916B (Pastor-Marazuela et al. 2021) 
and FRB20220319D (Ravi et al. 2023) have compa- 
rable upper-limits on local scattering and likely have 
DM, < 20pccm ?. The local (< 150 kpc) contribution 
to the DM of FRB 20220509G is less certain, due to the 
presence of both scattering and scintillation. We assume 
a flat prior between 0 and 100 pc cm" ?. 

For P(DMiem), we use the IllustrisTNG simulation 
to estimate the distribution of DM from the IGM for a 
source at a few hundred Mpc along similar sightlines to 
our FRBs. We take the functional form of DMiaw found 
by Zhang et al. (2021) for FRBs at z z 0.1, who also 
used IllustrisTNG. However, we independently estimate 
the mean DMijay for sightlines that intersect massive 
halos with Mgoo > 10135 Mo. We then exclude DM 
contribution from the halo itself to estimate the mean 
of DMizaw for that subset of sightlines. By doing this we 
account for correlations in the Universe's matter distri- 
bution: Sightlines that intersect clusters are more likely 
to intersect filaments and less likely to pass through 
voids, compared to typical positions. We find that for 
z S 0.1, 


(DMraw) x 22 pc cm? (xor) ‘ (8) 
where Da is angular diameter distance. This gives a 
mean value of 93pccm-^? with a 90% confidence inter- 
val on DMīcm = 48 — 260 pc cm? for FRB 20220914A. 
For FRB 20220509G we find a mean value of 76 pc cm? 
with a 9096 confidence interval on DMīıcm = 39 — 
175pccm-?. The estimated IGM DM distributions are 
shown in Figure 5 as dashed red curves. 

In these FRB directions, the Milky Way contribu- 
tion is taken to be normally distributed with mean 
70pccm^? and standard deviation 30pccm ?. This 
value includes the Galactic halo. These values are based 
on NE2001 (Cordes & Lazio 2002) and YMW (Yao et al. 
2017) and recent evidence that the Milky Way halo 
contribution to DM is smaller than previously thought 
(Bhardwaj et al. 2021; Ravi et al. 2023; Cook et al. 
2023). The two FRB sources are nearby on the sky, 
allowing us to use similar estimates for P(DMmw). 

We have attempted to make conservative estimates on 
the uncertainty of each component of the observed DM, 
which will result in a wider inferred DMjcy distribution. 
We warn that the various elements in the DM budget 
were estimated by different means and are subject to 
modelling uncertainties. With these caveats, the result- 
ing probability distribution of the ICM contribution to 
DM is shown in Figure 5. We find that Abell 2310's ICM 


adds 265-511 pc cm ^? at 90% confidence to the observed 
DM of FRB 202209144. In the frame of the cluster, this 
range is 295-568 pc cm" ?. These values make it unlikely 
that the host galaxy is near the front of the cluster, as 
such a high DMycm requires a significant path through 
the intracluster medium. In the case of Abell2311, we 
find a 9096 confidence interval of 16-172 pccm-^? for 
DMycm in the cluster frame. This range is consistent 
with an FRB embedded in a 1.6 x 1014 Mọ cluster at a 
projected offset of ~ 870 kpc. Still, it is difficult to say 
with certainty that the FRB's majority DM component 
is the ICM due to the presence of scattering and larger 
projected impact parameter of the host. 
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Figure 5. The relative probability distributions for each 
component of the observed FRB DMs. The distributions of 
DMicm are shown in the purple shaded region under the 
black curves. 


3.3. ICM magnetic field 


FRBs are the only known cosmological sources for 
which both DM and RM can be measured. If DMicm 
and RMjcm can be determined, FRBs allow us to mea- 
sure the mean line-of-sight magnetic field strength of 
a galaxy cluster, as shown in Eq. 6. The observed 
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RM of FRB 20220509G is —111.5 + 1radm"?, com- 
pared with the expected Milky Way foreground RM of 
—9+16radm~? in that direction (Hutschenreuter et al. 
2022). If the extragalactic RM is dominated by the ICM 
then the cluster-frame RM is ~ 120 rad m~?, and we can 
use P(DMiocw) to infer (By tom) at the impact param- 
eter of FRB 20220509G. We find that the average mag- 
netic field strength is 0.7-7.5 uG. Again, this is subject 
to the assumption that the observed RM is dominated 
by the ICM. 

The values we obtain for both RM and (Bj 1c) are 
in line with previous observations, as well as magneto- 
hydrodynamical (MHD) simulations of galaxy clusters 
(Carilli & Taylor 2002; Marinacci et al. 2018). Bohringer 
et al. (2016) analyzed the RMs from 1383 of extragalac- 
tic polarized sources and found values between —200 and 
+200rad m7? were common at Aog. The correspond- 
ing magnetic field strengths they deduced were ~ a few 
microGauss (Bóhringer et al. 2016). In Section 5 we dis- 
cuss how cluster FRBs from future surveys could better 
constrain magnetic fields in and around the ICM. 


3.4. Mean temperature from DM 


The ICM observables that we discuss in this work 
are DM, RM, X-ray brightness, and the SZ-y param- 
eter. Each is a weighted integral of n, along the line of 
sight. In Faraday rotation, the gas density is weighted 
by magnetic field strength, B;;. For X-ray luminosity, 
emission is the sum of n2 weighted by A(T;). In the 
case of the SZ-y parameter, ne is weighted by Te itself. 
X-ray derived density and SZ-derived pressure are com- 
monly used to estimate temperature (Eckert et al. 2013). 
Alternatively, the ratio of X-ray luminosity to ysz has 
been taken to constrain the gas temperature. For exam- 
ple, Planck Collaboration et al. (2016c) find the ratio of 
X-ray luminosity to ysz for the Virgo cluster at angular 
position 0 to be (ne A(T.)/(kpg Te)) (0). This quantity 
offers information on the cluster temperature profile un- 
der certain assumptions about ne(r) and A. 

We can attempt something similar using the ratio of 
ysz to DMrcw, which will be independent of ne. Eval- 
uating Eq. 2 and Eq. 3 at impact parameter b, we get, 


kp Te 
ysz(bi) ~ — 5 OT (ne) Lic (9) 
me C 
DMrom(b1) © (ne) Lrcm, (10) 


where Lrcy is the line-of-sight lengthscale of the cluster 
at that b1. By taking their ratio and re-arranging we 
can make a crude approximation for the mean plasma 
temperature at that physical offset, 


(11) 


=i 
E 7 YSZ DM 
= 1.15 x 107K (25) (am) (12) 


We estimate the projected y,, at each cluster’s b, from 
our fits to the Planck data described in Section 2. This 
gives ysz = (3.31 + 1.36) x 1079 for Abell2310 at 
520 kpc and ysz = (6.80 + 3.03) x 1077 for Abell 2311 
at 870kpc. Plugging these values into Eq. 12, we find 
(kp T.) m 0.8 — 3.9 keV for Abell2310, in agreement 
with the expected electron temperature for a cluster of 
similar mass at that impact parameter (Vikhlinin et al. 
2005; Pratt et al. 2007). This is the first time halo gas 
temperature has been measured using fast radio bursts. 
For Abell 2311 at 870 kpc we find (kp T.) ~ 0.6— 16 keV 
along that line of sight. The large range in tempera- 
ture is primarily due to the uncertainty in ICM DM for 
FRB 20220509G. 

While it is hard to match the statistical precision of 
an X-ray spectroscopic temperature measurement, there 
are benefits to using FRB DMs to derive gas tempera- 
ture: With FRBs we do not need to account for “clump- 
ing" effects due to the n2 dependence (Eckert et al. 
2015), the temperature will be approximately mass- 
weighted rather than the more complicated weighting 
intrinsic to spectroscopic measurements (Mazzotta et al. 
2004), and current X-ray facilities are generally not sen- 
sitive beyond ~ Rso9. The obvious drawbacks of using 
FRBs are that we have only a single sightline and the 
DM contribution of the IGM, Milky Way, and the host 
galaxy must be modelled. The former will be allevi- 
ated by FRB surveys with high areal density, providing 
multiple sightlines through individual halos. The latter 
may be aided by cross-correlation (e.g. (Madhavacheril 
et al. 2019; chi 2021)), which are less sensitive to host 
DM contamination. In Section 5 we describe how future 
surveys will allow for statistical FRB/cluster studies. 


4, ILLUSTRISTNG CLUSTER SIMULATIONS 


We use the IllustrisTNG simulation TNG300-1 (Nel- 
son et al. 2019) to estimate typical DM values of the 
ICM as a function of impact parameter. We seek also to 
estimate the variance between galaxy clusters of similar 
mass, and variance within a cluster for different position 
angle at fixed impact parameter. 

We first extract cylinders of length 300 Mpc and ra- 
dius 5 Mpc from the simulation box that contain massive 
halos. These allow us to estimate not only the ICM con- 
tribution to DM, but also the IGM along that sightline. 
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Figure 6. Four galaxy clusters from the TNG300-1 simulation whose masses are similar to those of Abell 2310 and Abell 2311. 
'The top row show gas surface density within 2 Mpc of the cluster center. The green and blue circles represent the impact 
parameter of FRB 20220509G and FRB 20220914A, respectively. The bottom row shows DMicm as a function of impact 
parameter, b, , with error bars for the two observed cluster FRBs. The solid purple line is the median DM over fifty azimuthal 
angles and the shaded region shows the 90% limits at that b,. The curves assume the FRB is embedded halfway in the ICM, 


ie. lpos = 0. 


We then calculate gas properties for each cell in that 
cylinder and estimate free electron density as, 


ne = fe pg XH /Mp, (13) 


where fe is the free electron abundance, p, is gas density, 
Xg is hydrogen abundance, and m, is the proton mass. 
We take Xj to be 0.75. The it cell's contribution to 
the total DM is, 


DM; = ni Ar; (1 + &) G (14) 


where Ar; is the distance through the cell at an impact 
parameter b;. The cell’s redshift is z;. For cellsize, Ax, 
Ar; = 24/ Az? — b?. The total DM is then XDM; for 
all cells that are intersected by the FRB. 

We compute DM profiles for all halos in TNG300-1 
that have Msoo > 10135 Mo. In Figure 6 we show four 
galaxy clusters with masses comparable to the masses 
of host clusters of FRB 20220509G and FRB 20220914A. 
The blue and green data points in the bottom row cor- 
respond to the estimated ICM contribution to the two 
FRB DMs, which agree well with the expected range at 
their respective impact parameters. There is significant 
variance within a cluster even at a fixed impact parame- 
ter, particularly beyond the virial radius (corresponding 
to © 1500 kpc given Rso9 ~% 800 kpc for these clusters). 


There is also some scatter between clusters, but the DM 
curves are roughly consistent. The mean values in DM 
vs. b, (solid purple lines in Figure 6) are also in agree- 
ment with analytic models for the cluster gas density 
(Vikhlinin et al. 2005; Prochaska & Zheng 2019). 

Compared to simulations of the CGM or gas in galaxy 
groups, cluster simulations are in relatively good agree- 
ment with one another as well as with observations (Op- 
penheimer et al. 2021). As the sophistication of cosmo- 
logical MHD simulations progresses, they will be an in- 
valuable tool for understanding the properties of FRBs 
that are impacted by cluster gas. This is especially true 
for magnetic field inference from FRB RMs and for gas 
temperatures outside of Rso9, as outlined in the follow- 
ing section. 


5. DISCUSSION 
5.1. Beyond the virial radius 


Outside of the virial radii of galaxy clusters, plasma 
is difficult to observe via X-ray. FRB DMs, however, 
are more sensitive to this gas due to the ne (rather than 
n2) dependence, and can constrain to ICM and WHIM 
on the outskirts of clusters. Fujita et al. (2017) have 
shown that by combining FRB DMs with the pressure 
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profiles from SZ observations, the temperature profile of 
the WHIM can be obtained beyond 1.5 Rago. 

In future FRB surveys, large numbers of localized 
FRBs will produce DM maps with high areal density, 
intersecting individual clusters along multiple sightlines. 
For example, the proposed DSA-2000 survey will spend 
roughly 596 of its time on deep-drilling fields that overlap 
with Rubin, XMM-LSS, Extended Chandra Deep Field- 
South, COSMOS deep fields (Hallinan et al. 2019b). 
This will lead to roughly 500 FRBs in just ~ 10 deg?. 
Nearly every massive galaxy cluster out to z — 0.3 will 
be intersected by at least one FRB sightline. We will 
have a wealth of multiwavelength data in this field, in- 
cluding SZ data and deep X-ray observations. Such a 
dataset will allow us to constrain the magnetic field in 
the ICM and the density profile in the WHIM, as well 
as gas temperatures outside of Roo. However, from 
Figure 6 we see substantial azimuthal and inter-cluster 
variance in DM at large radii. Analyses of the gas be- 
yond R599 ought to be done in tandem with cosmological 
MHD simulations. 


5.2. Prevalence of cluster FRBs 


Roughly 5% of sightlines from FRBs at z ~ 0.5 are 
expected to intersect galaxy clusters in the foreground 
(Prochaska & Zheng 2019). However, if an FRB sight- 
line has a projected offset from a cluster that is less 
than Roo, it is more likely that the FRB came from the 
cluster itself than from behind it, assuming the rate of 
FRBs is related to stellar mass and not just star forma- 
tion. This is because the overdensity of galaxies within 
massive halos is enough to counter the larger volume 
behind the galaxy cluster. Using R200 ^ 1 Mpc, the co- 
moving volume in region of a cone behind a redshift 0.25 
cluster out to z = 0.5 is ~7 x 10? Mpc?. Assuming a 
galaxy number density of roughly 107? Mpc^? (White 
et al. 2011), the volume behind the cluster should have 
only O(70) galaxies, compared with 10?-? galaxies in 
the cluster. We note again that this rough estimate as- 
sumes cluster galaxies are as likely as non-cluster galax- 
ies to produce FRBs. We have also ignored luminosity 
function considerations, which make FRBs behind the 
cluster to be more difficult to detect because they are 
less bright. It is therefore unsurprising that the first two 
localized FRBs that are known to be impacted by the 
ICM have come from member galaxies of clusters and 
not from behind clusters. 

We have cross-matched twelve sources localized by the 
DSA-110 with cluster catalogs and found two FRBs that 
reside in galaxy clusters. This suggests the fraction of 
FRBs from galaxy clusters is feFRB = DIT. using a 
9096 Poissonian confidence interval. We do not include 


other localized FRBs because we are not aware of con- 
certed efforts to cross-match those FRB positions with 
clusters. 

The large value of f.pgp has interesting implications 
for FRBs detected at other surveys, and on our interpre- 
tation of observables such as the FRB DM distribution. 
CHIME/FRB detects O(10°) FRBs per year, without 
sufficient localization precision to identify a host galaxy. 
'They may therefore already have detected hundreds of 
sources that reside in galaxy clusters. The host ICM 
could then dominate the high-DM tail of the total ob- 
served DM distribution (James et al. 2022), which would 
otherwise be interpreted as the most distant sources. 
Even with arcminute localizations, a careful statistical 
cross-match of CHIME/FRB sources with galaxy cluster 
catalogs could reveal this signal. 


5.3. FRB progenitor implications 


FRB 20220509G is the first source to belong to an 
early-type galaxy. This fits an emerging picture that 
FRBs can be produced in a variety of environments, 
including dense star forming regions, globular clusters, 
and pristine environments with little HII (Tendulkar 
et al. 2017; Marcote et al. 2020; Pastor-Marazuela et al. 
2021; Kirsten et al. 2022; Niu et al. 2022b). The full 
implications of FRB 20220509G's quiescent host will be 
fleshed out by Sharma et al. (in prep.). If more FRBs 
are localized to cluster galaxies with low rates of star for- 
mation, then progenitor models must explain evolution- 
ary channels that can produce FRBs in the absence of re- 
cent core-collapse supernovae. This would be promising 
for FRB applications research, as early-type galaxies are 
not expected to have significant magnetoionic plasma in 
their ISM or HII regions near the source, alleviating the 
problem of disentangling host DM, RM, and scattering 
from cosmological contributions. 


6. CONCLUSION 


We have discovered two fast radio bursts sources that 
reside in massive galaxy clusters. The host galaxy 
of FRB 20220914A is a member of cluster Abell 2310, 
whose ICM dominates the DM budget of the FRB. 
We have combined the DM of FRB 20220914A with 
SZ observations of Abell2310 to make the first esti- 
mate of the mean line-of-sight temperature of gas us- 
ing an FRB. FRB 20220509G belongs to an early-type 
galaxy (Sharma et al, in prep.) at a projected offset 
of 870 kpc from the center of cluster Abell2311. That 
cluster's ICM is found to contribute 16-172 pc cm? to 
the FRB DM, which is consistent with analytic mod- 
els and simulations of the ICM for a cluster with mass 
Mso0 © 1.5 — 2 x 101^ M. Polarization analysis of the 
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burst found significant Faraday rotation. Assuming this 
RM originates in the ICM, we constrain the mean line- 
of-sight magnetic field strength to be 0.75-7 uG in the 
intracluster gas. Roughly 1796 of our first sample of lo- 
calized FRBs were found to reside in galaxy clusters. 
While we expect this fraction to come down with time, 
galaxy clusters will likely play a significant role in up- 
coming FRB surveys, particularly in attempts at map- 
ping out the Universe's baryons. 
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